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Abstract 

We study the Dvali-Gabadadze-Porrati model by the method of the boundary effective 
action. The truncation of this action to the bending mode tt consistently describes physics in 
a wide range of regimes both at the classical and at the quantum level. The Vainshtein effect, 
which restores agreement with precise tests of general relativity, follows straightforwardly. We 
give a simple and general proof of stability, i.e. absence of ghosts in the fluctuations, valid for 
most of the relevant cases, like for instance the spherical source in asymptotically flat space. 
However we confirm that around certain interesting self-accelerating cosmological solutions 
there is a ghost. We consider the issue of quantum corrections. Around flat space 7r becomes 
strongly coupled below a macroscopic length of 1000 km, thus impairing the predictivity of 
the model. Indeed the tower of higher dimensional operators which is expected by a generic 
UV completion of the model limits predictivity at even larger length scales. We outline a 
non-generic but consistent choice of counterterms for which this disaster does not happen 
and for which the model remains calculable and successful in all the astrophysical situations 
of interest. By this choice, the extrinsic curvature K^ u acts roughly like a dilaton field 
controlling the strength of the interaction and the cut-off scale at each space-time point. At 
the surface of Earth the cutoff is ~ 1 cm but it is unlikely that the associated quantum effects 
be observable in table top experiments. 



1 Introduction 

In recent years several attempts have been made to formulate a theoretically consistent ex- 
tension of General Relativity that modifies gravity at cosmological distances while remaining 
compatible with observations at shorter distances. Well known examples in the literature 
are models of massive gravity ^ |2] , the GRS model [3| , the DGP model [H , and the new- 
born model of 'ghost condensation' These attempts are motivated both by the appealing 
theoretical challenge they represent, and, perhaps more interestingly, by the experimental 
indication that something strange is indeed happening at very large scales: the expansion of 
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the Universe is accelerating. This fact can be explained, as it is usually done, by invoking a 
small cosmological constant, or some sort of 'dark energy'. Nevertheless, it is also worth con- 
sidering the possibility that at cosmological distances gravity itself is different from what we 
experience at much smaller scales, in the solar system. This perspective could perhaps offer 
a new direction for addressing the cosmological constant problem (see for instance ref. [5]). 

In the present paper we concentrate on the DGP model. In this model there exists 
a critical length scale -Ldgp below which gravity looks four dimensional, while for larger 
distances it weakens and becomes five dimensional. To avoid conflict with observations, the 
parameters of the model can be adjusted to make Ldgp of the order of the present Hubble 
horizon Hq 1 . One cosmological solution of DGP is an ordinary Friedmann- Robertson- Walker 
Universe that, even in the absence of a cosmological constant, gets accelerated at late times by 
the weakening of gravity JTHH]. Moreover, even though the length scale of gravity modification 
is cosmological, one still expects small but non-negligible effects on solar system astrophysical 
measurements [HI El- This makes the model extremely interesting from the phenomenological 
point of view. 

In spite of these appealing features, the model is plagued by some consistency problems, 
as pointed out in ref. There it was shown that the theory has strong interactions at a 
tiny energy scale 

A ~ (Mp/L^Qp) 1 / 3 , (1) 

which corresponds to wavelengths of about 1000 km (a qualitatively similar, but quantita- 
tively different, result was found in ref. J2]). This means that below 1000 km the model 
loses predictivity, and one cannot perform sensible computations without knowing the UV 
completion of the theory. In the presence of curvature the situation seems even worse. In- 
deed in a non-renormalizable theory like the one at hand, quantum corrections necessarily 
imply the presence of an infinite tower of higher dimensional operators suppressed by inverse 
powers of the interaction scale A. A 'generic' choice of this tower of operators disrupts the 
classical computations of the original DGP model making it practically uncalculable, 
whenever the Riemann curvature is bigger than 0(1/Lq GP ). In practice the unknown UV 
completion of the model dominates physics above the tiny curvature 1 / L^ GP . In the case of 
a spherical source in asymptotically flat space this happens below a distance from the source 
given by the so called Vainshtein length |131 114j 

R V ~(R S L 2 DGP )^ , (2) 

where Rs is the Schwarzschild radius of the source. For instance, for the Sun Ry is about 
10 20 cm! Apparently one cannot compute the gravitational potential of the Sun at distances 
shorter than 10 20 cm. A completely equivalent problem has been pointed out in the case of 
massive gravity ^Hj. It should be stressed that, in the purely classical DGP model, Ry is 
precisely the length below which classical non-linearities in the field equations eliminate the 
phenomenologically unwanted effects of the extra polarizations of the graviton |13| I14j. The 
latter are usually referred to as the van Dam-Veltman-Zakharov discontinuity ^]. These 
non-linearities are therefore a virtue of the classical theory, but upon quantization they 
suggest the presence of many more uncontrollable interactions. 

Finally, the model possesses negative energy classical solutions, which signal the presence 
of instabilities already at the classical level. These negative energy solutions have however 
a typical curvature larger than l/L^Qp, so that their true presence, or absence, is fully 
dependent on the UV completion. 
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In this paper we try to understand to what extent all these results really threaten the 
consistency and the predictivity of the model. In spite of the above difficulties we will be 
able to draw fairly optimistic conclusions. In practice we will show that one can consistently 
assume a UV completion (= tower of higher dimensional operators) where the theory can 
be extrapolated down to distances significantly shorter than 1/A. Though our choice is 
consistent, in the sense that it is stable under RG evolution, it may at first sight look 
unnatural from the effective field theory view point. We will try to argue that perhaps it 
isn't. 

In the DGP model our world is the 4D boundary of an infinite 5D spacetime. One can 
integrate out the bulk degrees of freedom and find an "effective" action for the 4D fields. This 
was explicitly done in ref. where it was found that from the 4D point of view, beside the 
ordinary graviton, there is an extra scalar degree of freedom tt that plays a crucial role. This 
is essentially a brane bending mode contributing to the extrinsic curvature of the boundary 
like K^ v oc d^dyTi. It is this scalar that interacts strongly at momenta of order A. Indeed 
there exists a limit in which the strong interaction scale is kept fixed but all other degrees 
of freedom (namely the graviton and a vector N^) decouple. This means that in order to 
grasp the interesting physics related to the strong interaction we can restrict to the tt sector, 
which represents a consistent truncation of the theory. We will see that all the interesting 
features (good or bad) of the DGP model can be traced to the dynamics of this scalar tt. For 
our purposes we can therefore forget about the 5D geometric setup of the model and study 
a scalar theory with a specific cubic interaction, namely -^(dir) 2 Dtt. This largely simplifies 
the analysis. 

We start in sect. |3] by studying classical field configurations in this theory. For any given 
distribution of matter sources, the field equation for tt turns out to be a quadratic algebraic 
equation for the tensor of its second derivatives dJd v TT. Like ordinary quadratic equations, 
it possesses doublets of solutions: for any solution tt + (x), there exists a corresponding 'con- 
jugate' solution tt~(x) / tt + (x). By studying the stability of a generic solution against small 
fluctuations, one discovers that this conjugation exactly reverses the sign of the kinetic action 
of the fluctuations. This means that if tt + {x) is a stable solution, its conjugate tt~{x) must 
be unstable. It is therefore not surprising that classical instabilities were found in ref. 

Does the existence of these unstable solutions make the model inconsistent at the classical 
level? We argue that this is not the case: the classical stability of the tt system is guaranteed 
for a well defined class of energy-momentum sources and boundary conditions at spatial 
infinity, as we prove explicitly in sect.HJ In particular, we show that for any configuration of 
sources with positive energy density decaying at spatial infinity, and small pressure (smaller 
than p/3), the solution which is trivial at spatial infinity is stable everywhere 1 . Of course the 
conjugate configuration, which instead will diverge at infinity, will be unstable. The point 
is that the sets of stable and unstable solutions are disconnected in field space, they cannot 
'communicate'. It suffices that a solution is locally stable at one space-time point to ensure 
its global stability in the whole Universe, because there exists no continuous path of solutions 
that connects the stable and the unstable regions in field space. Unstable solutions simply 
correspond to the 'wrong' choice of boundary conditions. 

Unfortunately nothing prevents some interesting solutions from having the wrong bound- 
ary behaviour: it is the case of the de Sitter solution, which describes the self-accelerated 

If the sources have relativistic velocities an additional hypothesis is required, see sect. 0J 
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Universe we mentioned above. In ref. [IT] it was pointed out that such solution is plagued by 
ghost-like instabilities. In sect. El we check this result in our formalism, indeed finding that 
the self-interacting dynamics of tt is responsible for the existence of the de Sitter solution. 
The latter turns out to be nothing but the conjugate of the trivial configuration tt(x) = 0, 
and as such must be unstable. 

In sect.EJwe study the solution for tt in presence of a spherical heavy source. We reproduce 
in the tt language the results of refs. ^3E3G3E3) namely that the correction induced by the 
DGP model scales like (r/i?y) 3//2 relative to Newton's law. As stressed in ref. [S], the sign 
of this correction depends on the cosmological phase, i.e. on the behaviour of the Universe 
at infinity. In accordance with the general proof given in sect. |1J we explicitly check that 
fluctuations around this solution are stable. Also in this case the interesting features of DGP 
are encoded in the non trivial dynamics of the field tt, in a simple and transparent way. 

The results of sects. HJ El imply that for any astrophysical source and for a wide and 
relevant class of cosmological solutions no classical instability develops in the tt field. Finally, 
in sect. we discuss the problem of quantum corrections. We have already mentioned that 
a generic choice of higher dimension operators leads to uncalculability when the Riemann 
curvature is bigger than C(l/Lp GP ) ~ Hq. We look at this problem in more detail. We 
characterize higher order corrections, distinguishing between a classical and a quantum ex- 
pansion parameter. We notice that the large effects of the counterterms which were pointed 
out is ref. jllj are more related to the classical expansion that to the quantum one. We 
clarify this, by discussing first the structure of the 1-loop effective action. We show that 
the terms which are truly needed for unitarity and consistency, and which are associated to 
logarithmic divergences, do not grow in the regions of large curvature. In a sense this follows 
simply from the fact that the classical theory is well defined, i.e. there are no ghosts. Inspired 
by this simple example we then proceed to give a general and consistent characterization of 
the class of counterterm Lagrangians which do not destabilize the naive classical results. Our 
claim is based on a remarkable property of the tree level action for tt. The large curvature 
backgrounds around which one seemingly looses calculability correspond to a classical back- 
ground Id^d^TTd/ A 3 1 3> 1. Around such a background the quantum fluctuation ip = tt — tt c i 
receives a large positive contribution to its kinetic term 

£ ki „ = z^dPtpWtp) , (3) 

where is a non singular matrix proportional to d^d v TT c i/ 1\? itself. The interaction La- 
grangian remains however the same: no large factors there. The large kinetic term then sup- 
presses interactions (as seen for instance by going to canonical fields). At each spacetime point 
x one can define (suppressing spacetime indices) a naive local scale A(x) = ky/ Z(x) S> A, 
which is a function of the background field tt c [ itself. Our basic result is then the following. If 
the counterterm action depends on tt and A only via the running scale A and its derivatives, 
then the classical DGP action can be used down to the distance 1/A(x) at each point x. This 
implies, in particular, that the Vainshtein solutions within the solar system are always fine. 
On the other hand at the Earth surface one finds A ~ 1 cm -1 . In principle this could be 
worrisome, or interesting. In practice it is probably neither. This is because at the surface 
of Earth tt is well decoupled from matter. 

The counterterm structure we select must correspond to some specific UV completion for 
the physics at the 4D boundary. At the scale A we expect new states to become relevant. 
This scale is determined by the background expectation value of the massless field tt, or 
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more precisely d^dyir oc Ka U - Therefore the extrinsic curvature acts like a sort of dilaton 
controlling the scale of the new physics. A direct consequence of this property is that our 
ansatz can be motivated by an approximate scale invariance of the UV completion. 
We conclude with an outlook on possible developments of our results. 



2 Boundary effective action 



The DGP model |1] describes gravity in a five-dimensional spacetime M, with boundary dM. 
The action is postulated to be purely Einstein-Hilbert in the bulk, plus a four dimensional 
Einstein-Hilbert term localized at the boundary, with different Planck scales, 

S DGP = 2M| f d 5 xV^GR(G) + 2M| I d 4 x^R(g) - 4M 5 3 / d 4 x^K(g) , (4) 
JM JdM JdM 

where Gmn is the 5D metric, g^ u is the 4D induced metric on the boundary, and we have 
added the Gibbons-Hawking term at the boundary in order to obtain the 5D Einstein equa- 
tions upon variation of the bulk action |2()| . Matter fields and a possible cosmological constant 
are supposed to be localized on the boundary. A special role is played by the length scale 
-^dgp = l/ m = M4/M5: below -Ldgp gravity looks four-dimensional, while at larger length 
scales it enters in the five-dimensional regime. In order for the model to be viable, -Ldgp 
must be huge, at least of the order of the present Hubble horizon. 

In the following sections we will refer to the 4D boundary effective action obtained in 
ref. by integrating out the bulk degrees of freedom. We summarize here only the main 
results. By using spacetime coordinates {x^, y}, with the boundary sitting at y = 0, the bulk 
piece of the action can be rewritten in terms of ADM-like variables as 



'bulk 



oc 



2M| / d 4 x / dy V=g N [R(g) - K^K^ + K 2 ] , (5) 







where N = 1/yGvv is the lapse, = G yfl is the shift, g^ v = G^ v is the metric on surfaces 
of constant y, and the extrinsic curvature tensor is 

Kp, = ^ (d y g, w - V^N U - V„iV M ) . (6) 

One then expands the 5D metric and all the other geometric quantities around a flat back- 
ground, Gmn = f]MN + hMN, and integrate out the bulk in order to obtain an effective action 
for the 4D dimensional fields living on the boundary. The final result, at the quadratic level, 
is 

S bd y - Mj j d A x [\ h'^ u Uh'^ - \ ti Uti - m N'^AN'^ + 3m 2 vr Dtt] , (7) 

where A is a non-local differential operator, A = \f— □ = sj— r\^ v d^d u , and the kinetic terms 
have been diagonalized by defining 

h yy = -2A-K , N[ L = iV M - <9 M 7r , h'^ = h^ v - mir r]^ v . (8) 

By taking into account bulk interactions with higher powers of Hmn, one finds that the 
leading boundary interaction term is cubic in it, and involves four derivatives, 

= 'Ml [ d'x (dvr) 2 Dvr . (9) 
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The comparison of this interaction with the kinetic term of the tt field in eq. ((7J) immediately 
shows that the tt sector of the theory becomes strongly interacting at a very small energy 
scale A = (m?M 4 ) 1 /^ = M5/M4, which corresponds to a length scale of about 1000 km. 

(3) 

The boundary effective action we will use throughout the paper is the sum 5bdy + AS^ . 

(3) 

Notice that AS^ dy is the leading boundary interaction term in a quantum sense: it gives the 
largest amplitude in low energy scattering processes, and the scale A associated to it is the 
lowest of all strong interaction scales associated to further interaction terms. These terms 
are schematically of the form 

A£ bdy ~ M| d (N'T (dny (h'^y ~ mMl d [~^\ ( J^) ' , (10) 

where N'^ and tt are canonically normalized, and p + q + s > 3. Eq. Q corresponds to the 
term with p = s = and q = 3. 

However, in this paper we study classical field configurations, and from a classical point 
of view it is not immediately manifest why eq. © should be the most important interaction 
term, and why we should be allowed to neglect all the others. In fact, the strong interaction 
scale A has no meaning at the classical level (see discussion in sect. [7J. Nevertheless, we 
will check explicitly in the two specific cases we will deal with (the de Sitter solution and 
the spherically symmetric solution) that indeed all other interaction terms are sub dominant 
as long as the flat-space approximation is valid, i.e. as long as \h! \ <C 1. This fact can be 
understood in a general although formal way. First, if <C 1, in eq. ()10j) we can stick to 
terms with s = 0. Then, suppose we have a compact source for our fields: the region in which 
space is nearly flat is the region well outside the Schwarzschild radius Rs- We can mimic to 
be in this region by formally sending Rs to zero, i.e. by decoupling four-dimensional gravity 
from the source. In this limit the solution for the metric is trivial, h'„ u = 0. However, we 
do not want to end up in a completely trivial configuration, so we want to preserve the self- 
coupling of the tt field, as well as its interaction with the source: in terms of the canonically 
normalized tt field, the cubic self coupling is unchanged if we keep A = M 2 jM 4 fixed; in 
order to take into account the interaction with matter sources, we notice that the interaction 
Lagrangian between the 4D metric perturbation hav and the matter stress energy tensor 
is, by definition, ^h IJU/ T f " / . From the definition of tt, eq. (JSJ), we see that n interacts with 
matter via a term jMl^' Therefore, if we take the formal limit 

M 2 T 

M 4 , M 5 , Tfj, u -> 00 , — = const , -j^ = const , (11) 

we decouple 4D gravity while keeping the full Lagrangian for tt fixed. By 'full' we mean the 
sum of kinetic, cubic and source terms. We stress that when applying this limit to a spherical 
source of mass M, the Vainshtein radius Ry = RsLf^ GP = MM^/Mf remains fixed, showing 
that we are focusing on the genuine non-linearities of the DGP model. The point now is that 
in the above limit all coefficients of further interactions, eq. ()10|) . vanish, 

mMl d (-^kA " (3^ = ( T^-T 9 (NinOiry - . (12) 



m 1 / 2 M 4 \mM 4 J \M 2 J M 3 P /2+<?- 
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So far we neglected the dynamics of , but the above result tells us that we were justified 
in doing so: N'^ does not interact directly with matter, and in the limit we are considering all 
its interactions vanish, thus decoupling it from all other degrees of freedom. In conclusion, 
we see that in the limit in which the metric perturbation can be neglected, \h!^ v \ <C 1, eq. © 
gives the dominant interaction term also for classical field configurations. Conversely our 
approximation breaks down when \h! \ ~ 1, like when approaching a black-hole horizon. In 
fact when this happens we also have (dir/mM^) = dir ~ 1 (see sect. |HJ). Since tt roughly 
measures the bending of the brane, this simply tells us that our approximation breaks down 
when the brane can no longer be treated as approximately flat in 5D. In the final section we 
will briefly explain what we expect to happen in this regime. 



3 General properties of classical solutions 

We now move to study classical solutions in presence of matter sources on the boundary. We 
restrict to the tt sector, with cubic self-interactions and coupled to matter. As we showed in 
the previous section, this is a consistent truncation of the theory, in the sense that there exist 
a limit, eq. in which all other degrees of freedom decouple, and all further interactions 
vanish. 

From the results reported in the previous section we have that the full action for the tt 
field in flat space is 



S= (fx 



(13) 



-3(a*)>-^(0*)>D* + ^*T 

where T = T^, and we canonically normalized tt. 

The first variation of S with respect to tt gives the field equation 

6 Dtt - -L □(dvr) 2 + ^d,(8^ + Jl = , (14) 

which can be rewritten as 

3K->^ )2 + F (a)2 = -4^- (15) 

Notice that the terms involving third derivatives of tt have canceled out, thus leaving us 
with an algebraic equation for its second derivatives. In terms of the tensor field K liu {x) = 
— -^d^d^ir the above equation reads 

K 2 -(^) 2 -3K = - i ^-, (16) 

where K = K^n,. This result is easy to understand from a geometric point of view, by 
referring to the 5D setup of the model. The field equation that derives from varying the bulk 
action eq. © with respect to the lapse N is one of the Gauss-Codazzi equations, 

R{g) + K^K^ -K 2 = 0, (17) 
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which is an algebraic relation between intrinsic and extrinsic curvature of the boundary. 
The intrinsic curvature R(g) can be eliminated by means of the Einstein equations on the 
boundary, 

4M| (iV - \ 9ilv R) - 4M| [K^ - 9tMV K) = , (18) 

thus obtaining an algebraic relation between the extrinsic curvature K^ u and the matter 
stress-energy tensor. The extrinsic curvature is easily related to tt by expanding eq. ©, 



K, 



m 



Mi 

M 5 ' ~W 



dir 



+ q ( J_ . MiQh' 



+ o 



M 



(19) 

In the limit we are considering tt gives the dominant contribution, ~ — -^d^d u Tt = mK^ u , 
and this leads precisely to our eq. (fl6|) . This remark also tells us that eq. (fl6|) is more 
fundamental than its counterpart written in terms of tt, eq. (|15[). The latter is valid only 
in the limit of nearly flat space, while the former is a combination of a geometric identity 
and of the exact Einstein equations on the boundary, and as such is exact also in presence 
of large curvatures, as long as the classical theory makes sense. Of course in such a case the 
tensor K^ u is the extrinsic curvature of the boundary in units of m, rather than the second 
derivative of tt. 

Once a solution of eq. (jl(i|) for a given distribution of sources has been found, one may 
be interested in studying its stability. In order to do that, it is necessary to perturb tt and 
expand the action up to quadratic order in the perturbation <p, 



J (fx -3(dip) 2 - ^ (Qtt (dc/?) 2 + 2 d^TT d^y Utp) 



(20) 



By using the identity d^ipDip = d v [d v \pd fl ip — ^rj^ (d(p) 2 ], and integrating by parts, we find 



(fx 



d A x 



3(d<pY-2lKr,-ri ia ,K) d»<p&'<p 



(21) 



(22) 



This result deserves some comments. First, notice that the starting cubic interaction term 
involves four derivatives, but the resulting quadratic action for the fluctuation ip does not 
contain terms with more than two derivatives acting on the 99's. This remarkable feature 
makes the issue of stability more tractable, since it prevents the appearance of ghost states at 
high momenta, when higher derivatives contributions to the quadratic action could become 
important and destabilize the system. Second, the coefficients in S v depend only on the 
tensor K^ u , in an algebraic way. The fact that both the field equations and the condition 
for stability depend algebraically on K^ v will play a crucial role in what follows. 

A deeper insight into the structure of classical solutions comes from the following remark. 
For a given distribution of sources, eq. (jTK|) is a second-order algebraic equation for K^. 
Instead of having a single solution, it is reasonable to expect doublets of solutions K^, like 
in ordinary quadratic equations. Moreover, one expects that in general the stability of the 
system depends on the solution chosen inside this doublet, since the quadratic Lagrangian 
for the fluctuation tp, eq. (|22|1. depends on K^ u . Even better, given the simple dependence 
of the overall coefficient of d^ipd^ip in eq. (|22j) on K^ u , 



Z, 



-37/_ 



[IV 



2{K^ ~ V^K) , 



(23) 
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one can straightforwardly invert this relation and express K^ v in terms of Z^, 

K[iv = ^jiv ~ \{Z iiV — ^rj^uZ) , (24) 

where Z = Z^a- This is remarkable: now we can write our field equation, eq. (|16p. directly 
as an algebraic equation for the kinetic coefficient of the fluctuations, 

I Z *-( V f = 6 -W (25> 

In the above equation the linear term in Z^ has canceled out, thus leading to a purely 
quadratic equation. This explicitly shows that if Z+ u is a solution, so is its 'conjugate' 
7— = _ 7+ 

By definition eq. ()22j) now simply reads, 

S v = Jd*xZ ia ,(x)& l (pff'(p, (26) 

which shows that a solution is locally stable against small fluctuations if and only if the matrix 
Z^v, once diagonalized, has the usual 'healthy' signature (+,—,—,—). (In the following 
section we will be more precise in defining the requirements that ensure the stability of the 
system.) Therefore, if a solution is locally stable, its conjugate Z~ is locally unstable. 
The kinetic terms of the small fluctuations around the two solutions are exactly opposite. 

Coming back to K^ u , we conclude that if K^ v is a solution for a given matter distribution, 
then there always exists another solution K~ v = r]^ u — K^ u for the same matter distribution, 
and the stability properties of the two solutions are opposite. Of course, since by definition 
K^y = —jjd^d^-k, we have to check that K~ is indeed the second derivative of a scalar field. 
However this is automatic, since rj^ u itself is a second derivative, rj^ u = ^d iJL d u {x a x a ). 



4 Classical stability 

We now address in detail the issue of classical stability. We will show that under very general 
conditions a solution which becomes trivial at spatial infinity is stable everywhere. Before 
proving this fact we want to specify what we mean by 'stable'. In particular, we are going to 
characterize the stability of a generic classical solution by studying the tensor structure of the 
quadratic action of the fluctuations, eq. IJ2fi|) . around that solution. Notice that while for a 
stationary solution the concept of stability is perfectly defined, for a solution which evolves in 
time such a definition is a bit tricky. After all, if a solution depends on time the energy of the 
fluctuations around that solution is not conserved. For instance, one can expect to observe 
resonance phenomena associated to modes whose proper frequency is of the order of the 
inverse time scale on which the background solution evolves. A resonant growing fluctuation 
can be certainly seen as an instability, but its presence is not simply encoded in the tensor 
structure of the quadratic action. Our intention can therefore appear to be naive. However, 
the point is that our characterization is certainly sensible for ultraviolet modes, well above 
all possible resonances, on time and length scales much shorter than the typical time and 
length scales on which the background varies. In such a regime it is perfectly acceptable to 
study the stability of the system as if Zuv{x) were constant in space and time. We call this 
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concept of stability 'local', since it refers to the local structure of the quadratic Lagrangian 
of the fluctuations, and in principle nothing prevents the system from being stable in a given 
spacetime region (i.e. stable against short wavelength fluctuations localized in that region), 
and unstable in a different region. In this section we prove that under proper although 
general conditions this cannot happen. 

In a system described by a quadratic Lagrangian of the type of eq. 1)26(1 two physically 
different kinds of instabilities must be considered. The first is similar to a tachyonic insta- 
bility, and, roughly speaking, is related to the relative signs of the terms involving time and 
space derivatives. If one of the relative signs is wrong, there exists an exponentially growing 
mode. More precisely, consider the equation of motion deriving from eq. 1)26)1 . 

V d»d u ip = , (27) 

where, with the above discussion in mind, we neglected the space-time dependence of Z^ u . 
For a mode with four-momentum fc^ = (cj, k), this gives the dispersion relation Z tlu k^k u = 0, 
which can be cast in the form 

(Z 00 w + Z 0l k l f = (Z 0l Z 0j - ZwZn)k % W . (28) 

In order for the system to be stable, the solution for u> must be real for any wave vector 
k. The requirement is therefore that the 3x3 matrix Z^Zoj — Zq^Z^ is positive definite. 
Although not manifestly, this condition is clearly Lorentz invariant, since it corresponds to 
the requirement that the Lorentz invariant equation Z fJiU k^k u = has only real solutions 
(apart from an arbitrary overall phase multiplying k 11 ). 

The second kind of instability one has to consider is the ghost-like instability, which is 
related to the overall sign of the kinetic Lagrangian. Suppose that all relative signs are 
correct, so that all frequencies are real and there exists no exponentially growing mode. 
However, if the overall sign of the quadratic Lagrangian is wrong and our scalar field is 
coupled to ordinary matter, then it is possible with zero total energy to amplify fluctuations 
in both sectors, and this process is going to happen spontaneously already at the classical 
level. At the quantum level the problem is probably even worse, although it is not yet clear 
how to do sensible and consistent computations in a quantum theory with ghosts in the 
physical spectrum (see for instance |21l I22| b This kind of instability, unlike the previous, is 
there only in presence of interactions with other sectors. If the system is free of tachyonic 
instabilities, i.e. if ZoiZoj — ZooZij is positive definite, then the ghost-like instability is avoided 
if Z 00 > 0. 

In the Hamiltonian formalism both types of instability, although physically different, are 
formally on an equal footing. Their absence is guaranteed whenever the Hamiltonian is 
positive definite. For our action, eq. 1)261) . the Hamiltonian density of the fluctuations reads 

n(n, ip) = -n=-(n + 2Z 0t v i¥ >) 2 - Zij Vi^Vjip , (29) 

4Zoo 

where n is the conjugate momentum of (p, that is n = 2{Zqq ip — Zoi'Viip). 7i is positive 
definite for positive Zoo and negative definite Zij. The latter is a stronger requirement than 
the one we found above, since Z^Z^j is always a positive definite matrix. However here 
there is a subtlety. What must be positive definite is not the Hamiltonian density, but rather 
its integral in d 3 x, i.e. the total Hamiltonian. The mixed term ~ HViip that comes from 
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expanding the square in eq. (|29|) is conserved on the equation of motion, once integrated 
in d 3 x, since its time derivative turns out to be a total gradient. Therefore the integral 
H = J'd 3 xH , with 

w„(n, tp) = [|n 2 + (z 0i z 0j - ZooZ^v^Vjif] , (30) 

^00 

is conserved on the equations of motion of H = f d ?J xTL. Although Hq is not the Hamiltonian, 
a sufficient condition for stability (compact orbits) is that Hq be positive definite. Notice that 
this condition is weaker than that set by positivity of H, but it coincides with the condition 
found from the field equations. In summary, in order for a solution to be stable one needs 
that Zqq be positive and that the matrix Zq^Zoj — Zij be positive definite. 

Consider now a generic distribution of matter sources for the ir field. We stick to the 
case in which matter can be described by a fluid. Its energy-momentum tensor is thus of 
the form T^ IU = (p + p) u^u u + pr]^, where p is the rest-frame energy density, p is the rest- 
frame pressure, and it M is the fluid four- velocity. The tt field is sensitive only to the trace 
T = T^a = — (p — 3p), which does not depend on u^. We consider the case in which the 
sources are localized, i.e. their energy-momentum tensor decays at spatial infinity. In such a 
case far away from the source the trivial configuration = 0, Z^ u = —3p^ u is a solution of 
the field equations. Moreover, we start with the assumption that p = and p > 0, although 
we will relax it below. We will show that under these hypotheses the system is stable. 

Fix a generic x^ 1 . Z^ u is a (x-dependent) symmetric tensor, therefore we have some hope 
of diagonalizing it with a (x-dependent) Lorentz transformation. This is not always possible, 
since the Minkowski metric rj^ is not positive definite, and the spectral theorem does not 
apply. In the Appendix we show that Z^ u can be certainly diagonalized with a Lorentz 
transformation if the typical velocities of the matter sources are non-relativistic. This is 
plausible from a physical point of view, since, roughly speaking, the failure of the spectral 
theorem can be traced to a possible strong mixing between space and time that can lead to a 
zero- norm eigenvector of Z^ v . For relativistic motions of the sources we are not guaranteed 
that the diagonalization is possible. In order to proceed we have to assume it, so the result 
we are going to find can be applied only if such a diagonalization is possible. Let's therefore 
postulate that there exists a Lorentz frame in which Z^ v = diag(zo, z%, z 2 , z 3 ). In such a 
frame eq. reads 

F{z) = -\ [{zl + ... + zl)- (z oZl + z z 2 + ... + z 2 z 3 )] = 6 + > 6 . (31) 

In order for the system to be locally stable in x^ we want a positive Zqq and a positive definite 
ZfoZoj — Zij. In the Lorentz frame we are using this simply reduces to the requirement that 
all z^'s are negative. In the space of the eigenvalues of Z^ v we can associate to this condition 
four critical hyperplanes, namely z^ = 0, for p, = 0, . . . , 3. On these hyperplanes the system 
is marginally stable. We want to show that these critical hyperplanes do not intersect the 
set of solutions of eq. (|31|) . This result would permit us to conclude that the system is 
stable. In fact, at x — > oo the system is locally stable, since the K^ v = configuration is 
stable. It could become unstable in moving to finite x, but by continuity its trajectory in 
the eigenvalues space should cross one of the four critical hyperplanes. If the intersection 
between the critical hyperplanes and the set of solutions is null, this cannot happen. 
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The proof of the above fact is straightforward. Take a critical hyperplane, for instance 
zq = 0, and impose it as a constraint on the l.h.s. of eq. pTj) . One gets 



F(z)\ zo=0 = -§ [(z 2 + zl +zf) - (ziz 2 + Z1Z3 + z 2 z 3 )] 

= -i[( Zl -z 2 ) 2 + (z 1 -z 3 ) 2 + (z 2 -z 3 ) 2 ]<0, (32) 

while in order to solve eq. (|31j) F{z) should be larger or equal than 6. In conclusion, the z^s 
that lie on a critical hyperplane cannot satisfy the field equation. 

Of course, this argument does not prove that unstable solutions do not exist. Actually, 
as we showed in the previous section, for every locally stable solution there exists a locally 
unstable conjugate solution. What this line of reasoning proves is that the sets of locally 
stable and locally unstable solutions are topologically disconnected in the eigenvalues space, 
and therefore a solution which is locally stable at one point (infinity, in our case) is stable 
everywhere. And vice-versa: a solution which is locally unstable at one point is unstable 
everywhere. Notice that in the above proof the requirement that the sources are localized 
was needed only to ensure that at infinity the trivial configuration Kn V = is a solution 
of the field equation, so that there exists at least one point at which the system is locally 
stable. In the case of non-localized sources, if a specific solution is known by other means to 
be locally stable at some point, then its stability throughout the Universe is guaranteed by 
the same arguments as above. 

We have shown that, for sources with zero pressure and positive energy density decaying 
at spatial infinity, the solution which is trivial at infinity is stable everywhere. We can relax 
all these requirements, in order to find a condition under which the sets of locally stable and 
locally unstable solutions in the eigenvalues space come into contact, so that a solution could 
in principle be stable somewhere and unstable somewhere else. Allowing for generic p and 
p, eq. (|25|) now reads 

^>-+^f- < 33 > 

From this and eq. (|32|) we see that a critical surface intersects the set of solutions if and only 
if 

(p - 3p) < -12 A 3 M 4 = -12 m 2 M| . (34) 

This is a necessary condition for the development of local instabilities, starting from a solution 
locally stable at some point, for instance at infinity. 

So far we have dealt with very general features of the classical solutions of the DGP model. 
We now move to explicitly derive two interesting solutions and study their properties. 



5 de Sitter solution 

In refs. [7J |H] it was shown that the large-distance modification of gravity implied by the 
DGP model can make a Friedmann-Robertson- Walker Universe accelerate even in absence 
of a cosmological constant. In particular, one finds that at late time, when all the energy 
density and pressure have redshifted away, the Universe approaches an accelerating phase 
with constant Hubble parameter H = m. In other words, a solution of the DGP model 
is a four-dimensional de Sitter space with curvature radius L = 1/m, embedded in a five- 
dimensional Minkowski bulk. The stability properties of this de Sitter solution were analyzed 
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in ref. where it was shown that around this configuration the tt field has a ghost-like 
kinetic term. 

In this section we show that the four-dimensional effective action of the ir field, in absence 
of matter and cosmological constant, can account for this de Sitter solution. Of course, since 
de Sitter space is curved, while we are working in flat space, we have to restrict to a small 
patch around x = and look at de Sitter space as a small deformation of Minkowski space. 
This is sensible for x <C L, where L is the de Sitter curvature radius. 

Consider eq. (|16|) in vacuum, that is set T = 0. We are looking for a maximally symmetric 
solution, therefore we consider the ansatz K^ v = jKr]^. We obtain 

K = 4, K^ = r)^, tt(x) = -±A 3 ^ . (35) 

In the language of sect. |3J this is the solution which is conjugate to the trivial configuration 
Kav = 0. Therefore, we know that it must be unstable. The quadratic action of the small 
fluctuations <p around this solution in fact is 

S v = Jd 4 x [-3{8^) 2 - 2 (k^ - Vtlv k^ d^pd"^ = Jd A x3(d^) 2 , (36) 

exactly reversed with respect to the 'canonical' one. This is precisely the result of ref. jllj 
about the kinetic term of the ir field around the de Sitter background. 

To check that this solution corresponds locally to de Sitter geometry we must take into 
account the mixing between ir and the metric perturbation h^ u . This is straightforward to do 
by noticing that in terms of the Weyl-transformed metric h'^ v = h^ v — -^ttt]^ the quadratic 
Lagrangian, eq. Q, is diagonal, i.e. there is no mixing between it and h! . The solution we 
are studying corresponds to h' = 0, so that the metric perturbation is 

hfj, u = jj- TTVtiu = -\m 2 Vims x aX a . (37) 



This describes locally a de Sitter space. In fact, by performing a gauge transformation with 
parameter = \m 2 x^x a x a , the metric becomes 

hp U = hftv + <9 M e„ + d u e^ = m 2 x^x u , (38) 

and this is exactly the induced metric near x = on the hyperboloid 

v^^ + y 2 = A' ( 39 ) 

which describes the embedding in five dimensional flat space of a four dimensional de Sitter 
space with curvature radius L = 1/m = Ldgp- 

We can check explicitly that in the regime of validity of the solution we have found all 
higher order interactions, eq. (|1U|). can be consistently neglected. Our solution has n ~ A 3 
N' = and h' = 0, so that eq. (|1UI) reduces to 

A£ bdy ~ mMj d(m x) q ~ m 2 M 4 2 (m x) q ~ l , (40) 

to be compared with the cubic self-coupling of tt, which on the solution is of order A 6 x 2 = 
m 2 M 2 (mx) 2 . We can trust our solution only for x <^ Ldgp = 1/m, since we are working 
in a nearly flat space. In this regime eq. (|4Uj) gives a negligible contribution to the action. 
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6 Spherically symmetric solution 



R v 



(-Rs-^dGp) 1 ^ 3 ^ -^DGP , 



(41) 



As pointed out in refs. |3 E3> one could in principle test the DGP model by solar system 
experiments, at distances much smaller than the DGP scale, which instead is of the order of 
the Hubble horizon. This is because the correction to the Newton potential of an isolated 
body becomes important at distances of the order of the so-called Vainshtein scale, 

M \ 1/3 
i 2 M 2 ) 

where M is the mass of the body and Rg its Schwarzschild radius. At smaller distances 
the relative importance of this correction is small, of order (r / Ry) 3 / 2 , but nevertheless non 
negligible. An interesting fact is that, as shown in ref. 0, unlike its magnitude, the sign of 
this correction depends on the cosmological phase, i.e. on the behaviour of the Universe at 
cosmological distances. 

In our formalism all the relevant features of the DGP model are encoded in the self- 
interacting dynamics of tt. We will derive an explicit solution for the field generated by a 
point-like source, and show that its properties agree with the above results. 

Consider a static point-like source of mass M, located at the origin: T = —M6(x). We 
look for a static spherically symmetric solution ir(r), where r is the radial coordinates. In 
such a case the field equations further simplifies; in particular eq. (|14|) becomes 



1 - - 2 -> 

6 Vvr - -rrV(V7r) 2 + —Vvr V 2 vr 
A 3 A 3 



V 



6E - ^V\E\ 2 + -^EV ■ E 
A 3 A 3 



V 



6E + r 



A 3 r 



M 
2M4 



5 3 (x) 



(42) 



where we defined E = Vir(r) = fE(r), and we used V • E = ^^{r 2 E). After integration 
over a sphere centered at the origin, this reduces to an algebraic equation for E(r), 

4 E 2 \ M 



4Tir z 6E + 



A 3 



2M 4 ' 



whose solutions are 



A 3 



±J9r 4 + ±R v r-3r 2 



(43) 



(44) 



where Ry is the Vainshtein scale, eq. (|41|). The solution for tt is obtained by integrating E 
along r, 7r±(r) = J E±(r)dr. 

At small distances from the source the two solutions reduce to 



E±{r < R v ) 



A 3 4 /2 



The acceleration induced on a test mass by tt is roughly j^E. The relative correction to the 
Newton force is thus of order 



"4^2^ r 1 / 2 
1 

3/2 



(45) 



F w E/M 4 /j^V 

^Newton ~ M/M 2 r 2 ~ \R V J 



(46) 
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as expected. Notice that the magnitude of the correction does not depend on the solution 
chosen, but its sign does. This fact matches to what we mentioned at the beginning of 
this section, namely that the sign of the correction should depend on the behaviour of the 
Universe at cosmological distances. Indeed, although the two solutions E± are similar to 
each other near the source, far away they behave very differently: E + decays as 1/V 2 , while 
EL blows up as r 2 . Since E + at spatial infinity reduces to the trivial configuration K^ u = 0, 
one is tempted to conclude that E- should approach the de Sitter solution found in the 
previous section, since the de Sitter and the trivial solutions are conjugate to each other. 
However, this is not the case, as one can easily check by noticing that E- does not reduce 
to a maximally symmetric configuration at infinity. Indeed, E + and E_ are not conjugate 
in the sense of sect. El but rather they are 'spatially' conjugate, in the following sense. E + is 
constant in time, therefore it produces a purely spatial K^ v , that is K^ = 0. By the same 
arguments of sect. El one easily concludes that another static solution of the field equations 
is Kij = ^5ij — Kf-. E_ corresponds precisely to this solution. 

Since we are considering a localized source, with positive energy density, we know from 
the results of sect. El that the solution trivial at infinity, E+, is stable everywhere against 
small fluctuations. We can check this fact explicitly. For tt = 7r(r), eq. (|21|) in spherical 
coordinates reads 



J d 4 Xyfg 



3 + -to V 2 vr 
A 3 



(ip 2 - g ij diipdjip) + — ViV.-Tr^V 1 dkipdiip 



(47) 



where gij = diag(l, r 2 , r 2 sin 2 6) is the metric tensor. The interesting non-zero components 
of the connection are T T de = —r and T r M = — r sin 2 9, so that 



V r V r vr = E' , V e V e vr = r E , V^V^vr = r sin 2 6 E 
while V 2 vr = E' + -E. In the end we have 



(48) 



d 4 x^/g 



2 / . 2E 



3 + 



4 E 

A 3 r 



2 ( , E 

3+ A^r + 7 



(duff 



(49) 



where we defined (dnip) 2 as the angular part of (V</?) 2 , 



9 • 2 



(50) 



By direct inspection of the solution E + it is straightforward to see that the factors enclosed 
in square brackets are everywhere positive, thus proving the stability of the solution. 

Similarly to the de Sitter case studied in the previous section, also on this solution higher 
order interactions of the type of eq. (|10|) give a negligible contribution to the action, even well 
inside the Vainshtein region, where the extrinsic curvature is getting larger and larger. To 
see this, it is sufficient to plug into eq. (jTU|) the behaviour of the solution for Rs <C r <C Rv , 
namely vr ~ A 3 Ry 2 /r 1 / 2 , ti ~ R s /r and N' = 0. The result 



is 



AC 



bdy 



r J 



g/2+s+l 



(51) 
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which clearly shows that as long asr> R$ the dominant contribution comes from the cubic 
self-coupling of ir, i.e. from the term with q = 3 and s = 0. Like in the de Sitter case, as 
long as the space is nearly flat the leading interaction term is the cubic self-coupling of ir, in 
agreement with the general argument given in sect. |2j 

We want to conclude with a comment on the behaviour of the solution in the presence 
of several sources. Of particular interest for the discussion of the following section is the 
behaviour of the effective kinetic term coefficient Z^ v . In the Vainshtein region Zuv, the 
extrinsic curvature and the Riemann tensor R^ vpa scale as 

Z ~ KL DGP ~ VRL DGP ~ x/i^L DGP /r 3 / 2 , (52) 

where here on we neglect the indices. Notice, in particular, that the extrinsic and Riemann 
curvatures are associated to the same length scale: K ~ y/~R. For the asymptotic linear field 
the relation was instead K ~ RLdqp <C y/R. Now, the relation Z ~ y/R~ Ldgp is valid also in 
more general situations. For instance we have checked it explicitly in the case of a spherical 
matter distribution not localized at one point. It is also easy to derive it for the interesting 
case of two (several) pointlike sources, for instance the Sun and the Earth. As M© 3> M© we 
can neglect the Earth contribution in the asymptotic field. M© becomes relevant only close 
enough to Earth. In this region we can approximately treat the Sun field as a background 
giving just rise to a wave function Z© ~ y/M@ Ldgp / (M^r^ 2 ) where r© is the Earth-Sun 
distance. By defining a canonical field 7r© = ^/ZqTT, the field equation close to Earth is 
just eq. (|42|) with renormalized values A — > A© = Ay/Z@, M4 — ► M^Z & . Then applying 
the previous results we have that the renormalized Vainshtein radius for the Earth field is 
R Q y = / (M^h? Zq) = r^M®/M Q . This is precisely the distance from Earth at which the 
Riemann curvature becomes dominated by the Earth field. By applying the results of this 
section, at shorter distances the full if wave function is just Z ~ Zq(R qv /r) 3 / 2 ~ V^-^DGP; 
where R is now dominated by the Earth field. Notice in particular that Rq v is larger than 
the Earth- Moon distance, so that the effects of the tt field on the Moon rotation EI] can 
be studied by neglecting the presence of the Sun in first approximation. This is precisely 
like in GR: the Sun's tidal force can be neglected when studying the Moon's orbit in lowest 
approximation . 



7 Quantum effects 

Our discussion was until now classical. In order to characterize quantum versus classical 
effects it is useful to keep h ^ 1 and write the Lagrangian eq. (|13ft in terms of the non- 
canonical field 7r and the the 4 and 5D Newton's constants G4 and G5, 

C = -%d*f + Ud.?d\ + %T (53) 

(throughout this section we will neglect 0(1) factors and be schematic with the contraction 
of derivative and tensor indices). At the classical level physical quantities do not depend 
on the reparametrizations tt — > air and also C — > tjC, with a, ij constants. For instance, 
in the case of a localized source T = MJ 3 (x), the only invariant quantity is precisely the 
Vainshtein scale Ry = MG\jG±. Indeed we can also define a classical expansion parameter, 
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Figure 1: The graphs contributing to the 1-loop effective action, eq. 



measuring the size of classical non linearities, by taking the ratio of the second and first 
terms in the Lagrangian a c i = G^d 2 ir /G4 ~ K. Notice that a c j is invariant under the 
above reparametrizations. In the case of the spherical solution it is just a function of r/Ry. 
At the quantum level the Lagrangian becomes a physical quantity, £ — > r/C is no longer 
a 'symmetry', so that there is an additional physical length defined by Lq = UG\jG\ or 
equivalently an energy A = H/Lq. For Ldgp of the order of the Hubble length we have 
Lq ~ 1000 km. This analysis can be compared to the case of ordinary GR in 4D: the 
analogue of the Vainshtein scale is the Schwarzschild radius, while the classical expansion 
parameter is simply the gravitational perturbation hn V = gu,u~Vuf Similarly, at the quantum 
level the Planck length is the analogue of Lq. 

Like we have defined a classical expansion parameter a c \ we can define a quantum ex- 
pansion parameter a g = {Lqdd). In the case of a spherical potential we have a q = {Lq/t) 2 . 
Then we expect that any physical quantity A will be expandable with respect to the value 
Aq it takes in the linearized classical limit as a double series 

±=^ = J2anm^a%, (54) 

n.m 

with n, m positive integers. The dependence on only integer powers of a q follows from 
Lorentz symmetry. As an example of a quantity of interest we could consider the action 2 . As 
shown in ref. and as evident from eq. (|22|). the loops of the quantum fluctuation field ip 
generate terms involving at least two derivatives on the external background, i.e. involving 
K^ v = —d^dpTr/A 3 and its derivatives. In particular the tree level cubic interaction is not 
renormalized. For instance, cutting off the UV divergences at the scale A, the structure of 
the 1-loop correction is 

r i-i°°P = ^[a m A 4 + 6 m 5 2 A 2 + a 4 (c m lnA + 7 m )] (^rT • (55) 

The graphs contributing to r 1_loop are depicted in fig. ^ m counts the number of external 
legs, and by I m we indicate the finite parts. Notice that r 1 " 10013 / ' C d = r 1 ~ loop /(5vr) 2 has 
the structure of eq. (|54jl. We should stress however that the power divergent terms are 
strictly speaking not calculable, in the sense that they fully depend on the UV completion. 
On the other hand the log divergent terms are also associated with infrared effects and 
are thus calculable. They correspond in eq. to powers of n that are multiple of 3: 
o? ~ ft 793 <9 6 = E e /A 6 is the genuine loop counting parameter. For instance by working in 



2 Notice that although the action is not a physical quantity at the classical level its relative variation 5S/Sq is 
physical, and indeed defines a c i. 
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dimensional regularization we would only get the minimal set of divergences that are required 
for the consistent definition of the theory, and these are precisely the logarithmic ones. 

Now the important remark is that when the classical parameter a c i is larger than 0(1) 
the expansion in eq. Q54|) breaks down even if a q <C 1. What basically happens is that large 
classical effects pump up the quantum corrections, and the result becomes fully dependent 
on the series of counterterms. Both power and log divergent terms at a c i > 1 formally swamp 
the classical action for a large enough number of external legs (to in eq. (|54|l). This implies 
that the details of the UV completion of DGP in general dominate the physics for a c / > 1 
or at distances from a point source that are shorter than the Vainshtein radius Ry. Then 
we should not even be able to compute the gravitational potential at solar system distances, 
although they are still much bigger than the genuine quantum length Lq. A similar problem 
has been pointed out in the case of massive gravity |15j . 

There are however indications that the situation may not be necessarily as bad as it 
seems. First, as we have already remarked, the problem arises when expanding in the num- 
ber of external legs, that is when expanding in the classical background field. For any given 
amplitude with fixed number of legs the loop expansion, sum over n in eq. (|54|). is convergent 
as long as r > Lq. Second the classical theory we are dealing with makes perfect sense even 
below the Vainshtein scale as we have proven in the previous sections. For a well defined class 
of energy-momentum sources and of boundary condition at spatial infinity, no pathology or 
ghostlike instability is encountered in the ir system. Indeed even if the interaction term 
involves four derivatives, the quadratic action around the background remarkably contains 
quadratic terms that have only up to two derivatives: this prevents the appearance of ghost- 
like modes. This should be compared to what happens when higher derivative corrections 
to a kinetic term become important and ghost states appear: in those cases it is mandatory 
to assume that extra quantum effects eliminate the problem, or better to declare the regime 
of validity of the theory to be limited. Based on these two remarks we want to consider the 
structure of the quantum effective action more closely. As we did previously we decompose 
7r into classical background and quantum fluctuation: tt = ir c i + ip. The tree level Lagrangian 
for the fluctuation is 

C v = Z^&'ip&'tp + ^d^cpd^Dcp (56) 

and Z^v, see eq. (|23|) . contains all the dependence on the classical background tt c [. In the 
regions where classical non-linearities are large, the eigenvalues of the kinetic matrix become 
large, Z ~ K 3> 1. Now, as long as Z does not lead to ghosts, a large Z will enhance the 
gradient energy and suppress quantum fluctuations in tp. This somewhat contradicts our 
previous conclusion that things are out of control when K S> 1. In order to simplify the 
discussion and the notation consider a toy model where the background matrix reduces to a 
scalar Z^ v = —r]^ u Z. As long as Z varies slowly over space, it makes sense to define a local 
strong interaction scale A(x) = A.\J Z(x), which describes the local scattering of quantum 
excitations p. Of course this scale is obtained by writing the action in terms of a locally 
canonical field tp = tp ■ y/Z, and the definition makes sense as long as the typical distance by 
which Z changes is itself much bigger than 1/A(x). It is easy to calculate the divergent part 
of the 1-loop effective action r 1_lo °P by working with the canonical field tp. The dependence 
from Z in r 1_loop comes from two sources: 1) the Jacobian of the rescaling in the path 
integral, 2) a spacetime dependent mass term —m 2 (x)tp 2 generated by the kinetic term after 
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4 V Z 2 / 2 Z 

The contribution from the mass term corresponds to the usual Coleman- Weinberg effective 
If W = ^ (AfWm 2 - (m 2 ) 2 In A UV + f ? W ) , (58) 



rescaling, 

2 i ( d ll zd»z \ id 2 z 

m ( x ) = 7 ( ^ 2 ) " • (57) 

ibution from the n 

potential 

I 

where the finite part w is expected to be a fairly complicated non-local expression that 
fully depends on the spacetime behaviour of m 2 . For instance, in a central field where Z oc r a , 
like in the Vainshtein region, we have m 2 ~ 1/r 2 and we expect r^ w = — (m 2 ) 2 lnr plus 
terms of order m 4 , but without log enhancement. On the other hand the Jacobian factor 
yields formally 

rJ = -l5 4 (0)lnZ. (59) 

This term does not depend on the derivatives of Z, and should correspond to a quartic 
divergence. We can make this explicit by regulating the delta function as 5 4 (0) = A^y/ (2tt) 4 . 
Notice that in principle the ultraviolet cutoff Auv can be as high as the local scale A(x). To 
comment this results let us focus first on the logarithmic divergent term, as it is completely 
unambiguous. We see that when Z > 1 this term does not blow up, as one could have 
concluded by looking at the series expression eq. . The effects of a large Z are resummed 
in the effective mass m 2 , in a form that depends only on the size of the relative derivatives 
d n Z/Z. For instance in the Vainshtein region the logarithmic part leads to a correction to 
the classical action which is of order 

' 1 ln(A uv r) . (60) 
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This should be compared with the classical action which is of order A 4 (Ai2y) 2 (iiy/r) 5 / 2 = 
r 2 A 6 (r): we find that the quantum correction is subdominant for r > r c = l/(A 4 Ry) = 
LqI \M&Rs)- Indeed r c is precisely defined by the equality r c = A _1 (r c ) as one should have 
expected. Notice that for the interesting cases of macroscopic sources r c is much shorter 
that Lq = 1/A = 1000 km. We will have more to say on this below. We have more 
comments to make on this simple 1-loop computation. Eq. (|58|) . when expanded in powers 
of K = Z — 1, corresponds to the calculation of the c m coefficients of eq. (|55|). This object 
describes the 1-loop RG evolution of the effective Lagrangian. On the other hand the initial 
conditions for this evolution at some scale, for instance A, are determined by arbitrary 
local counterterms contained in the I n . In general, as long as locality and the symmetries 
are preserved, there is little restriction on these terms from the effective Lagrangian view 
point. In particular they could become large just when Z is large. For instance a function 
of the form Z N (dZ/Z) 4 would be a perfectly acceptable counterterm, with a perfectly fine 
analytic and local expansion around the point K = Z — 1 = 0. For JV > 1 a term like 
this would dominate the action below the Vainshtein scale. What distinguishes a damaging 
term like this from the mild effects that we have explicitly computed? The difference is 
very simple: this term depends on the absolute magnitude of Z, while the others don't. 
After the rescaling, the dimensionful parameters in the tree level (p Lagrangian are given by 
A = k^fZ and by the derivatives dZ/Z = 25A/A, d 2 Z/Z = <9 2 A 2 /A 2 , etc. The operator in 
question cannot be written just in terms of these local scales since Z N = A 2N /A 2N and the 
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original scale A does not influence local quantities in the tree approximation. We can then 
characterize in a fairly simple way the class of 1-loop counterterms with four derivatives (log 
divergences) for which perturbativity is preserved below the Vainshtein scale (in fact below 
Lq): it is given by local polynomials of A, d 2 A, etc., weighted by the appropriate powers 
of A to match dimensionality. Nothing particularly interesting is learned by considering the 
structure of the power divergences. After all this is not surprising since they are totally 
scheme dependent. However the lesson we drew from the logarithmic ones is clear enough 
to allow us to formulate a simple consistent requirement, a conjecture, valid for all type of 
counterterms, and to all orders, such that the quantum corrections remain under control 
way below the Vainshtein scale. The requirement is that the only dimensionful quantity 
determining the local counterterms be the local scale A(x), 

£ct = A(x) 4 F (|) , (61) 

where F schematically represents any infinite polynomial in dA, d 2 A, etc. Notice that the 
above counterterm Lagrangian should be added with A written in terms of the physical field 
7T = 7r c i + (f, as this decomposition is arbitrary. Then we need to check that eq. 1)61)1 does 
not lead to new interactions for the quantum field ip associated to scales lower than A. If 
that were the case our ansatz would not be not self consistent since beyond 1-loop it would 
introduce a new scale in the effective action for the background field. However it is easy to 
verify that this consistency check is satisfied. Notice that the minimal set of counterterms 
that are needed in DR with minimal subtraction belongs to this class. However our definition 
is more general. For instance a term with no derivative A 4 is consistent with eq. (|61|) but 
is not generated in DR, since it corresponds to a quartic divergence. The physical meaning 
of our ansatz is also clear. The scale A truly represents the scale at which new degrees of 
freedom come in. The UV divergences in the low energy effective theory are cut off at this 
scale. For instance the vacuum energy is essentially A 4 up to small derivative corrections 
associated to the space-time variation of A. From the point of view of the fundamental UV 
theory the field Z = 1 + K is a sort of dilaton that controls the mass gap A. At the points 
where Z goes through zero the effective description breaks down since the gap vanishes and 
new massless states become relevant. This is another way to argue that the points where 
Z is negative, the ghost background, are outside the reach of our effective field theory and 
should be discarded. 

We based all our discussion on the simplified case Z^ u = —Zr)^, in which a rescaling of 
ip allowed us to define the running scale A(x). In the general case this cannot of course be 
done. The procedure then is to work at each point with a short distance expansion of the ip 
propagator, 

where the dimensionless coefficients (functions of x) di,... are calculated in perturbation 
theory around a slowly varying Z background. In spite of this more complicated tensor 
structure the discussion is qualitatively the same: the 1/Z in the propagator arrange with 
the 1/A 3 in the vertices to yield a sort of running scale A(x). However the dynamics at the 
scale A will not be Lorentz invariant. Anyway we stress that in the simple interesting case 
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of a spherical source inside the Vainshtein region we have 



V - (Rv/rf/ 2 A i 



(63) 



where A^ u is a constant matrix. Then a simple rescaling of the quantum field is sufficient to 
study the effective action. 

At this point one will wonder about naturalness. Are we advocating an incredibly tuned 
structure of counterterms not dictated by any symmetry but just by our desire to save the 
model? From the point of view of the effective theory expanded around tt = it looks like we 
are making a big tuning. This was the viewpoint of ref. for DGP, and, earlier, of ref. ^5] 
for massive gravity. However, as we have proven in the first part of this paper, the validity of 
the classical theory extends well beyond the point tt = 0. By validity of the classical theory 
here we include the fact that A never crosses 0. All we have done in this section is to demand 
uncalculable UV effects not to drastically perturb this picture throughout the patch in field 
space where the classical theory works. The fact that this patch is sensibly bigger than the 
point 7r = leads to stronger constraints on the counterterms. But when looking at the 
region tt ^ there is more to say. Around the points with ddn 3> A 3 our ansatz can indeed 
be motivated by a symmetry: scale invariance 3 . In the limit dd7r 3> A 3 we can neglect the 
quadratic kinetic term and work just with the cubic Lagrangian at tree level. The cubic 
Lagrangian is obviously scale invariant by assigning tt dimension zero. We can also trivially 
absorb the powers of A in the ir field: ir = Att. The background tt c i spontaneously breaks 
scale invariance. In particular it generates the mass scale A 2 ~ ddix c i. Then it is obvious that 
our ansatz preserves scale invariance. The original quadratic kinetic term now reads A 2 (9fr) 2 
representing a soft explicit breaking of scale invariance. When treating this term as a small 
perturbation it is then natural to generalize our eq. ()61|) to include all the terms proportional 
to positive powers of A 2 /A 2 . Now, one should be aware that scale invariance, unlike other 
global symmetries, cannot be readily used to naturally enforce parameter choices. Basically 
this is because conformal invariance in the interesting cases is never exact in the far UV. For 
instance quantum gravity breaks it. This limits the use of conformal symmetry as a substitute 
of supersymmetry in attacking the gauge hierarchy problem (see for instance a discussion in 
ref. [2S])- However one can imagine a situation where the scale M at which a models flows 
to a conformal point is much smaller than the physical and conformal breaking cut-off Atjv- 
Then even though at the fixed point there may exist relevant deformations their coefficient 
will be naturally small (in units of M) if the corresponding operator is irrelevant at the 
UV scale. An example could be a 4-fermion term suppressed by l/A^y in the microscopic 
theory, which becomes relevant at the fixed point thanks to a big anomalous dimensions. 
Such an example suggests that a theory may remain naturally scale invariant in a wide range 
of scales below the (conformal breaking) cut-off. In our case conformal symmetry is surely 
broken at the scale M5, but we do not need to go that high in energy. This keeps open 
the possibility that our ansatz may indeed be technically natural. It would be interesting to 
investigate these issues in more detail, perhaps trying to come up with an explicit realization 
of a confer mally invariant completion. 

Finally we want to discuss the physical implications of our results. The running scale 
close to a spherical source can also be written as 




3 



(64) 



3 We thank Sergei Dubovsky for waking us up on this. 
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from which the critical length r c at which the classical background ir c i stops making sense is 
just r c = 1/RsMq. We have r c < R$ for R$ > I/M5 ~ 1 fm, so that the quantum effects in 
the field it have no implication on any relevant astrophysical black hole. Of course we still 
have to study the classical effects of the field it on 4D black holes. 

In the presence of several macroscopic sources, like in the solar system, eq. (|64[) is dom- 
inated by the one generating the biggest gravitational tidal effect, i.e. the one dominating 
the Riemann tensor ~ i?s/r 3 . At the surface of the Earth it is the Earth field itself that 
dominates A and we find 

A" 1 ~ 1 cm . (65) 

In principle this is the length scale below which we cannot compute the gravitational potential 
at the surface of the Earth! Notice that this scale exceeds by almost 2 orders of magnitude 
the present experimental bound on modifications of gravity |241 125[ 126] . Notice also that 
we cannot play with M5 in order to increase A ffi by 10 2 . This is because A ffi oc \[M^ while 
Lpjgp ^ so that we would lower at the same time Ldgp down to the unacceptably 

low value of ~ 1 kpc. However the 1 cm scale we are discussing is the one at which tt, not 
the usual graviton, goes into a quantum fog. The graviton at this scale is still very weakly 
coupled. Moreover the field n at distances just a bit bigger that 1 cm is a very small source 
of macroscopic gravity, its contribution relative to the ordinary one being of order 

3 

'^V s £~10-«. (66) 
RvJ Ajjg 

Then it is reasonable to expect that whatever UV physics will take over below 1 cm it will 
not raise the effects of this sector on the gravitational potential by 16 orders of magnitude 
right away. Nonetheless this result may be viewed as yet another motivation to study the 
gravitational forces below a mm. 



8 Outlook 

We conclude by briefly mentioning some directions in which our analysis could be extended. 
First of all, there is the issue of stability of time-dependent classical solutions. As we have 
already stressed, our proof of classical stability is rigorously valid only for matter sources 
with non relativistic velocities. It would be interesting to understand to what extent the 
stability is guaranteed in the presence of relativistic sources. Perhaps one can encounter 
astrophysical or cosmological situations in which the kinetic Lagrangian of the fluctuation 
around the background passes through zero in some spacetime region. Such a region would be 
characterized by a local scale A = 0. This clearly indicates a breakdown of the effective theory, 
and in principle one expects unsuppressed quantum effects to take over, with interesting 
phenomenological signatures. Of course these are, by definition, out of the reach of the 
effective theory, and we can say nothing about them without a UV completion of the model. 

On the other hand it would also be interesting to study whether the approximate confor- 
mal symmetry of our ansatz can tell us more about the UV completion. Notice in passing 
that by neglecting the original quadratic kinetic term, the 7r-sector becomes a critical 4 ghost 

4 The kinetic term is not negative but zero. 
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condensate model. It is the presence of a classical energy-momentum source, and not the 
7r self-interaction, that causes the ir to condense, giving thus rise to a well defined effective 
field theory. 

Another point which may require extra study concerns 4D black holes. Are they different 
from those of ordinary GR? The fact that the field ir, when approaching the Schwarzschild 
radius, gives a more and more negligible correction to the Newton force suggests that black 
holes should be unaffected by its presence. However the very expansion we use breaks down 
close to the Schwarzschild radius, as shown by eqs. ()1(J1 151|) . Therefore a (numerical) study 
of the exact solution may be required. Of course the problem we are encountering here is 
just that close to the horizon we cannot consider the brane as approximately flat. A better 
viewpoint is gained by considering the full Einstein equations at the boundary in eq. (|18j) . 
In the region R$ <C r <C Ry the extrinsic curvature term, which controls deviations from 
GR, represents a small perturbation. This can be phrased as 

K 



R{9) -^^-L^ = W « V/W^ ~ 7? . (67) 

corresponding to the fact that the geometry is approximately Schwarzschild. The only way 
this state of things can drastically change is if the extrinsic curvature develops a singularity 
while approaching the Schwarzschild region. Without an explicit computation we cannot 
exclude this possibility, but we find it very unlikely. In analogy with GR, the horizon should 
not be a special locus from the point of view of the curvature. 

Finally, it would be worth analyzing in more detail the connection between DGP and 
massive gravity, and investigating the possibility that our results, or some of them, can be 
exported to the latter. This is far from trivial, since massive gravity apparently doesn't share 
with DGP two remarkable features that turned out to be crucial in our analysis: the fact 
that the equation of motion for the strongly interacting field depends algebraically on the 
coefficients of the quadratic action around the background, and the absence in that action 
of terms with more than two derivatives. 
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Note added 

In a recent paper ref. |23| it was also argued that the UV completion of DGP should not 
destroy predictivity. Such a claim is in the same spirit of our paper, although ref. [2Z1 is 
mostly based on an interesting u-model version of DGP, where there are just scalar fields. 
In that toy model there is no analogue of our tt field, with all its associated peculiarities, so 
that it is not immediately clear to us that there is a direct connection to our work. It would 
be interesting to further investigate this possibility. 
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A Diagonalization of in the non-relativistic limit 

In this Appendix we want to prove that whenever the matter sources have non-relativistic 
velocities the tensor Z^ v is diagonalizable with a Lorentz transformation. 

First, from the very definition of Z^, eq. (|23|) . we see that if is diagonalizable, so is 
Zav- Given Kav in a generic frame, by means of a spatial rotation we can always align the 
spatial vector Koi to the x axis and cast Kav in the form 



( 


ko 











hi 



















V 












(68) 



Now, the point is that if by means of a boost along x we succeed in diagonalizing the upper- 
left 2x2 block we are done, since then by a spatial rotation we can diagonalize the full matrix. 
The problem is thus reduced to the diagonalization of a 2x2 symmetric tensor with a boost, 
i.e. to find two real eigenvectors of the matrix 

*°'-(Z~*)- (69) 

orthonormal with respect to the 2D Minkowski metric r\ a g (now the indices a and (5 run over 
and 1). The eigenvalues of K a g are 



A± = l(k! - k ) ± l^Jih + k ) 2 - 4A& , (70) 

which are real and distinct as long as the off-diagonal element is small, \koi\ < ^\k\ + ko\- 
In such a case there exist two real eigenvectors, orthonormal with respect to rj a g. When 
instead \koi\ — ► g |&i + ko\ the two eigenvalues become equal, there is only one eigenvector 
and it has zero norm. Beyond this point, when |fcoi| > h\ki + ko\, both the eigenvalues and 
the eigenvectors are complex. 

Now we want to make contact with the physics. The tensor K^ u is related to the second 
derivatives of n, = —^d^d u Tr, so that 

ko = -jjir , k i = -jsd x 7r , ki = -jjd^TT . (71) 

If the matter sources are stationary we can certainly choose tt to be independent on time. 
In such a case ko = koi = 0. Now we can look at the case in which the sources are evolving 
with non-relativistic velocities as a perturbation of the stationary case. In particular, if the 
typical velocity of the sources is v, then we expect that the time derivatives of it will be 
suppressed by powers of v with respect to the spatial gradients, namely 

ko ~ v 2 k\ , fcoi ~ v k± . (72) 

In the non-relativistic case, » < 1, the two eigenvalues in eq. (|7U|) are certainly real and 
distinct, 

X + = kl -^ + 0(v 3 k 1 ) , A_ = -ko + ^ + 0(v 3 h) . (73) 

ki k\ 

This shows that for non-relativistic sources the 2x2 matrix K a g is diagonalizable with a 
Lorentz boost, so that the full 4x4 is diagonalizable, and so is Z^. 
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